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L0 INTRODUCTION 

The study of the burning of a single droplet is an ideal problem from which to gain fundamental understanding of diffusion 
flame characteristics. Droplet combustion is a complex physico-chemical process that involves a chemically-reacting two-phase flow 
with phase changes and yet simple experiments and analysis can be used to attain important insights into the burning rate, flame 
dynamics, kinetic extinction and disruption processes. It is a subject that has been actively studied for the past 40 years with most 
of the fundamental experiments being performed under reduced-gravity conditions for direct comparisons with 
theoretical/computational analyses that invoke spherical symmetry assumptions. 

In the earlier studies, the effects of sooting on the overall burning characteristics were not considered [Kumagai and Okajima, 
1957; Kumagai, Sakai and Okajima, 1971]. However, recent microgravity investigations performed at the NASA-LeRC droptowers 
(Droplet Combustkm Experiment of Professor F.A. Williams and Professor F.L. Dryer) and others indicate that effects of soot 
and sootcloud formation may be significant during the lifetime of the droplet and therefore must be included in the analysis [Knight 
and Williams, 1 980; Shaw et al., 1988; Choi et al., 1990; Jackson et al., 1991; Avedisian, 1994]. 

1. 1 Droplet Burning Rate 

Sooting mechanisms can have profound effects on the overall droplet burning rate. In an earlier study on microgravity n- 
heptane droplet combustion experiments, vaporization rates that were as much as 30% to 40% lower than the classically-accepted 
results were reported [Choi et al., 1990]. The main conclusion drawn from that investigation was that the effects of formation and 
accumulation of the sootshell and its asymmetric configuration caused by small degrees of relative droplet/gas convection were 
responsible for the observed discrepancies. There are several mechanisms through which soot/sootshell formation can affect the 
burning rate. Soot formation is an endothermic process and will act as a heat sink for the energy that would otherwise be used to 
vaporize the fuel. Although, exact figures for the conversion to soot is not presently known for microgravity droplet combustion, 
results from normal gravity pyrolysis experiments suggest relatively large conversion rates - 30% [Vranos and Liscinsky, 1984]. 
Furthermore, heat loss at the flame due to soot radiation can reduce the effective transfer number by decreasing the amount of energy 
available to vaporize the fuel. Related to these phenomena is the accumulation of the formed soot in the region bounded by the droplet 
and the flame. This represents unrealized chemical energy resulting in the reduction of the heat of combustion available from the 
vaporized fuel. 

1.2 Extinction 

Analysis of the extinction mechanism during droplet combustion is important for extracting overall chemical kinetic rates 
for large hydrocarbons for which detailed information does not exist [Chung and Law, 1 986; Williams and Dryer, 1 994] . Extinction 
occurs when the flow time of the reactants becomes smaller than the time needed for the reactants to be consumed (when the 
Dahmkohler number is reduced below a critical value). The critical Dahmkohler number is given by where d^ is the extinction 
diameter, D g is the diffusivity and t c is the characteristic chemical reaction time. 





(i) 


However, extinction can occur prematurely if radiative losses from the flame decreases the flame temperature and therefore affect the 
chemical reaction time. Chao and coworkers [Chao et al., 1 990] suggest that radiative losses from the flame can promote the onset 
of kinetic extinction. For droplet flames that produce soot, it is likely that soot radiation will dominate and therefore affect the kinetic 
extinction process [Grosshandler, 1992]. Thus far, only a limited number of extinction measurements have been reported for n- 
heptane droplets burning in air [Yang and Avedisian, 1988; Hara and Kumagai, 1990]. However, the extinction diameter is usually 
very small and not always reproducible under mildly sooting conditions. Without accurate determination of the soot concentration 
near the flame, it is not known to what degree soot formation may eventually affect droplet extinction. 
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1.3 Disruptive Burning 

In a previous work by Knight and Williams [1980], disruption was observed during the combustion of alkane droplets in 
a microgravity environment Disruption is a process by which the original droplet is violently shattered and burns as smaller satellite 
droplets. In subsequent experiments in the NASA-Lewis droptower [Shaw et al., 1988], it was found that the likelihood of observing 
the disruption behavior was higher for fuels (and cm- conditions such as higher oxygen concentrations) that produced high soot loading. 

Shaw and W illiams [1990] presented results that suggest that the absorption of heavy soot precursors on the droplet surface 
can create ranHitims in which an outer layer (of soot precursors) can reach temperatures exceeding the superheat temperature of the 
original fuel. Another possible mechanism suggests that the soot particles can come in contact with the droplet and cause rapid 
heating Although a definitive explanation does not exist, it is clear that the sooting mechanism is very important in the disruptive 
burning process. 

1.4 Thermophoresis 

For droplet combustion, soot particles form in the fuel-rich region of the diffusion flame through pyrolysis reactions. These 
reactions are high-activation energy processes and thus maximum soot production will likely occur near the flame [Shaw et al. 1988]. 
Once soot is formed in the fuel-rich region, they are acted upon by viscous drag (caused by Stefan flow) and Thermophoresis. 
Thermophoresis is a mechanism through which particle transport is caused by a temperature gradient in the surrounding gas-phase 
[Friedlander, 1977], Waldmann and Schmitt [1966] analyzed the viscous drag and thermophoretic force acting on a small spherical 
particle by calculating the momentum transferred per unit time to the particle by the impinging and rebounding gas molecules. At 
equilibrium [£F is equal to 0], these two forces are balanced. The resulting equation can then be solved for the flux, pV t 


„ -3 pdT/dr , , (2) 
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The term, pV„is defined as the thermophoretic flux. The thermophoretic flux is the convective flux necessary for the above-mentioned 
force equilibrium on the particle to be maintained. 

Opposing the thermophoretic forces is the viscous drag caused by the Stefan flux, where (dm/dt) is the droplet mass 
gasification rate. 
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Eqn. 3 is valid under conditions of negligible surface regression rate compared to the gas velocity at the surface which is satisfied when 
the ratio of the liquid to gas densities is large. The sootshell will form at the radial position where the viscous drag and thermophoretic 
forces are in balance. 

The above formulation is valid for soot particles that are much smaller than the mean free path. Jackson and coworkers 
[1991] implemented the thermophoretic model of Talbot and coworicers [1980]. In this model, the sootshell location can be calculated 
by assuming the soot particle dimension. Therefore, the measurement of the characteristic dimensions of the soot particles (including 
the radius of gyration and primary particle size using thermophoretic sampling and TEM) will be important in comparing the 
experimental and theoretical sootshell locations. 

1.5 Droplet Flame Dynamics 

The flame diameter is also affected by the soot formation through reductions in the specific volume [Choi et al., 1993]. After 
the initial period of fuel vapor accumulation (during which fuel vapor fills the increasing volume between the droplet surface and the 
flame front), quasi-steady conditions should exist under which the gasified mass will equal the mass consumed at the flame front. 
However, under sooting conditions, portions of the gasified fuel are converted to soot (which possess negligible specific volume when 
compared to that of the gaseous species). Therefore, soot accumulates in the sootshell and never reaches the flame front. 

2.0 The Need for Accurate Measurement of Soot Volume Fraction and Temperature 

Since the pioneering microgravity droplet combustion experiments of Kumagai and coworkers back in 1 957, there have been 
numerous theoretical, computational and experimental studies analyzing the burning characterstics of isolated droplets. In all 
subsequent microgravity droplet combustion studies, the degree of sooting was estimated by visual observations of the sootshell. This 
practice could not provide detailed information regarding the transient distribution of the soot and temperature in the region bounded 
by the flame and the droplet. For all of the important parameters of interest in droplet combustion including the variations in the 
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burning rate, promotion of extinction, reduction in flame size, onset of disruptive burning and soot dynamics, it is clear that accurate 
measurements of the degree of sooting and the temperature field within the region bounded by the droplet and the flame is needed to 
further advance our understanding. 

In the present study, the sooting behavior will be studied using light extinction measurements and Abel deconvolution to 
determine the radial distribution of soot volume fraction; two-wavelength optical pyrometry to determine the soot temperature within 
the region between the flame and the droplet; and thermophoretic sampling/transmission electron microscopy to determine the soot 
morphology (radius of gyration and primary particle size). In addition, computational modeling [using the ChemicaUy-Reacting-Flow- 
Model of Princeton] of the soot particle dynamics involving the balance between thermophoresis and Stefan drag will be compared 
with the experimental measurements. The modeling efforts will be advanced interactively with the experimental measurements of 
the soot volume fractions and the soot particle dimensions. 

These experiments will be performed for a wide range of conditions to vary the sooting propensity of the droplet using 
various fuels (including decane, heptane, toluene/methanol mixtures, benzene/methanol mixtures and ethanol etc.), pressure (0. 1 to 

1 .0 atm), oxygen indices (0.15 to 0.50), initial droplet diameters (1 .0 to 5.0 mm) and inerts (N 2 , He, C0 2 , etc.) as parameters. 

3.0 EXPERIMENTAL APPROACH 

3. 1 Soot Volume F raction Measurements 

Kadota and Hiroyasu [1984] first measured soot concentration in droplet flames using optical techniques. In their study, light 
extinction measurements were performed at a fixed location below the droplet and the radial variation was deduced by assuming that 
the regressing flame maintained constant spherical structure. This procedure may produce realistic information since the flame 
geometry is nearly spherical in the region below the droplet However, the same approach will produce unrealistic results in the region 
above the flame where it displays both axial and radial variations. 

In the present experiment, the soot volume fraction will be measured using a full-field light extinction technique that was 
developed at NASA-Lewis Research Center by Greenberg, Griffin and Ku [Griffin and Greenberg, 1 992; Greenberg and Ku, 
1 995] for eventual use in measuring the soot volume fraction for microgravity droplet combustion. Feasibility of this technique was 
investigated by performing soot volume fraction measurements for suspended droplets burning under normal-gravity conditions. 

Fig. 1 displays the schematic of the experimental apparatus similar to the system of Greenberg and Ku [1995]. Light from 
the HeNe laser is first passed through a variable neutral density filter and expanded through a spatial filter/collimator to produce a 
50mm diameter beam. The beam is directed through opposing optical ports (each fitted with 5 1 mm diameter quartz windows) of the 
stainless steel droplet combustion chamber. The transmitted beam is imaged using a high-resolution CCD camera (with attached 
interference and neutral density filters). The entire burning history is recorded onto a S-VHS recorder and subsequently the image 
is digitized frame-by-frame using a high-resolution (512x480 pixels) data acquisition board. Measurements indicate that the 
laser/detection system is very stable with temporal fluctuations of the light intensity less than 0.5% of the total magnitude of the 
intensity. 

For the normal-gravity droplet combustion experiment, a single droplet of 2mm diameter is deposited on a bead that is 
formed at the end of the suspending quartz fiber (diameter is approximately 150 micrometers) using a hypodermic needle. The 
orientation of the fiber as shown in inset A in fig. 1 was designed to minimiz e interactions with the burning by reducing its contact 
with the flame (For the proposed microgravity experiment, droplets will be suspended by using cross-hairs of 50 micron silica carbide 
fibers). The droplets are ignited using an induction coil discharge igniter. 

For each horizontal plane, the line of sight projection data, Pj/x) 9 is given by: 


P f (x) = ffj(x 2 +y 2 )dy 


(4) 


The Abel transform allows the calculation of the field distribution f y (r) corresponding to the projection data with the underlying 
assumption that it is radially symmetric [Santoro et al., 1983; Dasch, 1992; Sunderland et al., 1995]. The soot volume fraction 
distribution (based on extinction) is found from the measured transmission of each projection: 

T t (x) = y-(x) = exp[ ~ff v \ Z (* 2 + y 2 ) d y] (5) 

Aft A* 

lo -• 

where AT, is the dimensionless extinction constant of soot In this study, K e value of 8.8 which was measured for soot in the post-flame 
region of a rich acetylene flame was used [Choi et al., 1995]. Equation 5 can be rearranged into the form of the projection value as: 
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For a typical experiment, nearly 200 images are produced, therefore, full-field determination of the soot volume fraction measurements 
requires the analysis of more than 30,000 projection data files. To increase the efficiency, a full-field algorithm was written to analyze 
the entire experiment [Greenberg and Ku, 1 995 have also written a similar algorithm]. 

In the present study, the soot volume fraction is measured as a function of radial position for hydrocarbon mixture flames 
at various heights above the droplet. Figure 2 displays the resulting G/G 0 (where G corresponds to graylevel) profiles for t = 2 sec 
after ignition The entire image is filtered using a 3x3 low-pass filter and die individual vertical position data is 5-pt averaged. The 
horizontal axis corresponds to the column position and the vertical axis corresponds to row positions. The maximum attenuation is 
approximately equal to 0.4. From this profile, the soot volume fraction is determined using the full-field Abel-deconvolution algorithm 
[See fig. 3]. A curve-fitting algorithm is used to define the center of the sooting region for each vertical position and the soot volume 
fraction on either side of the center is averaged to obtain the final output. The soot volume fraction increases from as a function of 
height and radial position from the center. The soot volume fraction is maximum near the flame at all vertical locations which is 
consistent with the structure of diffusion flames. 

3.2 Soot Temperature Measurement Technique 

The temperature of the soot particles within the region between the droplet surface and the flame front will be measured using 
a two-wavelength pyrometry and Abel deconvolution [Faeth, 1 994]. The emission from the flame will be simultaneously directed 
to two separate CCD cameras (with a narrow bandpass filters at 800 nm and 900 nm) using a beamsplitter. Due to the small size of 
soot particles, the discrepancy with the local gas temperatures is expected to be small [KOylU and F aeth, 1993]. The emission intensity 
at each central wavelength measured for a vertical plane corresponds to the projection value of the integrated emission. 
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where I bx is Planck’s function. The projections values for 800 nm and 900 nm can be deconvolved separately using the Abel 
technique to determine the corresponding product of (K x I bx ) as a function of radial position. The ratio of the emission at 800 nm 
and 900 nm can then be used to determine the temperature using the two-wavelength pyrometry formula [Sivathanu and F aeth, 1 990] . 

3.3 Soot Morphology Measurement Technique 

The morphology of soot aggregates will be determined by performing thennophoretic sampling experiments [Megaridis and 
Dobbins, 1 988]. Thin ribbons of copper mesh grids will be rapidly inserted into the flame to obtain instantaneous distribution of soot 
morphology. Temporal variation of the soot will be determined by inserting different ribbons at pre-determined times. Subsequent 
analyses will involve obtaining transmission electron micrographs and analyzing the radius of gyration using a non-subjective digital 
image processing technique that was developed as part of this project [Choi and Zhou, 1 995]. 

4.0 PROGRESS TO DATE 

a) A droplet combustion chamber which utilizes several features of the DCE droptower hardware is being fabricated at NAS A- 
LeRC to take advantage of the design solutions. The UIC team will incorporate all of the optical diagnostic equipment 
to the experimental package. It is anticipated that an operational soot volume fraction/temperature rig to be used in the DC-9 
and in the 2.2 sec droptower will be ready by mid-summer of 1995. 

b) A full-field normal-gravity droplet combustion apparatus to perform soot volume fraction measurements was implemented 
at UIC . Performed experiments using various hydrocarbon mixture droplets. 

(results will be presented at the Central States Sectional Meeting of the Combustion Institute, 1 995). 

c) Developed an interactive Abel-deconvolution algorithm to determine soot volume fraction for full-field light extinction data. 
This program manipulates the full-field digitized image data and produces the I/Io information at every single pixel location. 
The information is then used to determine the soot volume fraction information for each horizontal plane. 


Similar program was also written for two-wavelength pyrometry. 

d) Developed a digital image processing program to determine the radius of gyration from transmission electron micrographs. 
The students are fully trained in TEM and Selected Area Electron Diffraction Techniques. 

e) Submitted a proposal (with P. Ferkul and D. Schultz of NASA-LeRC) to perform soot volume fraction measurements for 
droplets using the Glovebox facility. 

f) The PI was at Lewis during the summer of 1 994 working with Dr. Randy Vander Wal using laser- 
induced Incandescense for performing soot volume fraction measurements in droplet combustion. 

(results were presented by RVW at the Eastern States Sectional Meeting of the Combustion Institute, 1 994). 

5.0 FUTURE WORK 

a) Perform soot volume fraction/temperature measurements under normal-gravity for test conditions outlined in section 2. 

b) Incorporate optical diagnostic apparatus on the droplet combustion apparatus (that is being built at NASA-LeRC) for use 
in 2.2 sec droptower and aboard the DC-9. Reduced-gravity droplet combustion experiments for soot volume fraction and 
temperature measurements are scheduled for mid-summer, 1 995. 

c) Collaborate with Princeton researchers in incorporating revised thennophoretic models to calculate sootshell location. 

d) Determine soot extinction constants for various fuels to be used in soot volume fraction experiments. 
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